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Abstract 

We compare the structural and mechanical properties of static packings composed of fric- 
tionless convex (ellipses) and concave (rigid dimers) particles in two dimensions. We employ 
numerical simulations to generate static packings and measure the shear stress in response to 
applied simple shear strain as a function of the aspect ratio and amount of compression. We find 
that the behavior near jamming is significantly different for ellipses and dimers even though both 
shapes are roughly characterized by the aspect ratio and possess the same number of transla- 
tional and rotational degrees of freedom per particle. For example, we find that ellipse packings 
are hypostatic (not isostatic as found for dimers), display novel power-law scaling of the static 
linear shear modulus and contact number with the amount of compression, and possess stress- 
strain relations that are qualitatively different from that for dimers. Thus, we observe that 
important macroscopic properties of static packings of anisotropic particles can depend on the 
microscale geometrical features of individual particles. 

1 Introduction 

Significant progress has been made in understanding the jamming transition that occurs in collec- 
tions of frictionless spherical particles with purely repulsive short-range interactions, where systems 
transform from liquid- to solid-like states as a function of increasing packing fraction or decreas- 
ing applied shear stress [1, 2, 3]. Key findings include the power-law scaling of the static shear 
modulus with packing fraction above the onset of jamming [4, 5], the identification of a growing 
lengthscale as the system approaches the jamming transition [6], above which the system can be 
described as an elastic material [7] , and an abundance of low-energy excitations in the density of 
vibrational modes [8]. Much of this behavior stems from the fact that frictionless, static packings 
of spherical particles are typically isostatic since they possess the minimal number of contacts per 
particle ziso = 2d, where d is the spatial dimension, required for mechanical stability [9]. 

However, less progress has been made in understanding the jamming transition in particulate 
systems composed of nonspherical particles, despite the fact that these systems display striking 
mechanical [10] and rheological [11] properties, and are more relevant for industrial applications and 
in nature. An important difference between static packings of frictionless spherical versus ellipsoidal 
particles is that the latter are typically hypostatic, not isostatic, with fewer contacts than required 
to constrain all of the translational and rotational degrees of freedom using straightforward counting 
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Figure 1: Ensemble averaged (a) contact number zj and (b) packing fraction (pj at jamming as a 
function of aspect ratio a for dimers (squares) and ellipses (circles) for = 480 particles. 



arguments [12, 13]. Previous studies have found that for ellipse packings the contact number at 
jamming zj < Zjso = 2(ij, where = 3 in two dimensions {2D), over a wide range of aspect 
ratios a [12]. We include similar results from our simulations of static ellipse packings in Fig. 1 (a). 
zj < Ziso for small aspect ratios, but slowly approaches a value Zj that is only a few percent below 
the isostatic value as a increases. The packing fraction at jamming (pj for ellipses, shown in panel 
(b), possesses a peak, which is only a few percent lower than the crystalline value for spherical 
particles, near a ~ 1.5 [12, 14]. 

2 Motivation 

In this manuscript, we investigate the generality of these results for the behavior near jamming 
of frictionless, anisotropic particles by comparing the structural and mechanical properties of two 
classes of nonspherical shapes: convex (ellipses) and concave (rigid dimers) particles. We find 
that the behavior near jamming for rigid dimers differs significantly from that for ellipses even 
though both shapes are roughly characterized by the aspect ratio and possess the same number of 
translational and rotational degrees of freedom per particle. 

We find several key differences between the structural and mechanical properties of static pack- 
ings of dimers and ellipses. First, our simulations indicate that static packings of dimers are isostatic 
(not hypostatic as found for ellipses) with z ~ Z[so contacts per particle over the full range of aspect 
ratios studied as shown in Fig. 1 (a) ^. Second, ellipse packings display novel power-law scaling 
of the static linear shear modulus G and contact number z — zj with (j) — [15] — both scale 
linearly with (f) — cpj. In contrast, for dimer packings G and z — zj scale as {(f) — cpj)^'^, which is 
the same scaling found for static packings of spherical particles [4]. Third, we find that the shear 
stress-strain relations for packings of dimers and ellipses are qualitatively different. For example, at 
large compressions, the stress response (below the yield stress) to applied strain depends strongly 
on aspect ratio for ellipses, but it is nearly independent of aspect ratio for dimers. Also, at small 
compressions, dimer packings display nearly perfect plastic response in a region of strain where 
ellipse packings possess roughly linear response. 

This manuscript is organized as follows. In Sec. 3 we describe the computational methods 
for generating static packings of ellipses and dimers and then applying quasistatic simple shear 
to measure the mechanical response. In Sec. 4, we present our results for the linear static shear 

^It is still an open question which concave particle shapes yield isostatic packings, and which do not. 
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Figure 2: Definition of the aspect ratio a = a/b (ratio of the major to minor axes) for (a) ehipses 
and (b) dimers. 

modulus, contact number, stress-strain relation, and particle rearrangement statistics. In Sec. 5, 
we discuss our conclusions and identify possible future research directions. We also include five 
appendices, which provide the details necessary for calculating the packing fraction for dimers, 
contact distance between ellipses, and forces, torques, and stress tensor for anisotropic particles [16]. 

3 Computational methods 

We performed computational studies to measure the structural and mechanical properties of static 
packings of rigid dimers and ellipse-shaped particles in 2D. The particle shapes are shown in Fig. 2. 
The rigid dimers are formed by fusing identical disks together. We study aspect ratios a = a/b in 
the range 1 < a < 2, where a and b are the length of the major and minor axes, respectively. To 
inhibit crystallization, we focus on bidisperse mixtures of particles: 2N/3 particles with minor axis 
b and larger particles with minor axis 1.46. The particles are enclosed in square simulation cells 
with box length L and periodic boundary conditions. System sizes were varied from 24 < < 480. 

The particles interact via soft, pairwise, purely repulsive linear spring potentials. The total 
potential energy is therefore given by 



where fij is the vector separation between the centers of particles i and j, e is the characteristic 
energy scale of the interaction, Q{x) is the Heaviside function, and aij is the contact distance that 
in general depends on the orientation of particles i and j, fii and /ij, and fij. 

Contact distance Determining the interactions between dimers is straightforward: one can 
identify overlaps between individual disks (monomers) on different dimers. Thus, the contact 
distance between disk i on a given dimer and disk j on a different dimer is aij = (6j -|- bj)/2, and 
the total potential energy can be obtained by summing up the contributions V{rij/aij) over all 
disk-disk interactions for disks on distinct dimers. 

The contact distance (Tij between ellipses is more difficult to calculate than that for dimers. 
We define Uij as the distance at which two ellipses will first come into contact when moved along 
their center-to-center direction while their orientations are held fixed. Fig. 3 illustrates how aij 
is measured for ellipses i and j with orientations fii and fij at separation fij. We calculate the 
contact distances Gij analytically in systems of bidisperse ellipses using the Perram-Wertheim for- 
mulation [17, 18, 19, 20]. Further details are provided in Appendix B. 

Packing-generation algorithm We generate static, zero-pressure packings of bidisperse dimers 
and ellipses using a generalization of the compression/decompression method employed in our 
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Figure 3: (a) Ellipses i and j with orientations /i, and jlj at center-to-center separation fij. (b) 
The contact distance fjjj is obtained by identifying the point of contact when the two ellipses are 
brought together at fixed orientation. 

previous studies of spherical particles [21, 22]. Representative packings of dimers and ellipses are 
shown in Fig. 4. We briefly outline the packing-generation procedure here for completeness. 

We begin the packing-generation process by choosing random initial particle positions and 
orientations within the simulation cell at packing fraction (po = 0.20 (which is well below the 
minimum packing fraction at which frictionless packings of ellipses and dimers occur in 2D for 
1 < Q < 2). We successively increase or decrease the minor axes of the particles while maintaining 
the aspect ratio, with each compression or decompression step followed by conjugate gradient 
minimization [23] of the total energy (1). The system is decompressed when the total energy at 
a local minimum is nonzero — i.e., there are finite particle overlaps. If the potential energy of the 
system is zero and gaps exist between particles, the system is compressed. The increment by which 
the packing fraction (/> is changed at each compression or decompression step is gradually decreased. 
The process is stopped when the total potential energy per particle V/eN -C 1. Further details of 
the packing-generation algorithm are provided in Appendix C. 

The packing fraction cpj, contact number zj, and mechanical response are used to characterize 
each static packing. The packing-generation process is repeated at least 100 times at each a to 
generate configurational averages. Once the packings at jamming onset are generated, they can be 
successively compressed by small amounts A(p, followed by energy minimization at each step, to 
yield sets of configurations at fixed (j) — 

Quasistatic simple shear To determine the mechanical properties of static packings of dimers 
and ellipses, we studied their response to quasistatic simple shear at fixed area. We first initialized 
the system with an unstrained packing at a given (j) — <j)j and successively applied to each particle 
i small affine simple shear strain steps ^7 = 6x/L along the x-direction with a gradient in the 
2/-direction: 



where fi = {xi,yi) is the location of the center of mass of particle i. To be consistent with simple 
shear, at each strain step, the angle 9i = cos~^ {jli • x) that particle i makes with the x-axis was 
also rotated: 



Xi ^ Xi + 6jyi 



(2) 



Oi cot (cot 9i + 57). 



(3) 
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Figure 4: Snapshots of static packings of = 240 bidisperse (a) rigid dimers and (b) ellipses with 
a = 1.5. 

Each shear strain step was followed by conjugate gradient energy minimization using Lees-Edwards 
(shear periodic) boundary conditions [24]. For most studies, ^7 = 10~^ with accumulated strains 
7t = 1. We verified that smaller strain steps yielded similar results. During the quasistatic shear, we 
measured the shear stress (in units of e/6), contact number, and statistics of particle rearrangement 
events. These measurements are described in the following Sec. 4. Details of the stress calculations 
are provided in Appendices D and E. 

4 Results and discussion 

We present several measurements of the structural and mechanical properties of ellipse and dimer 
packings as a function of aspect ratio and compression (p — (pj including the contact number, shear 
modulus, yield stress, and other features of the stress-strain relations. 

Contact number at jamming The contact number is defined hy z = Nc/{N — Nr), where 
Nc is the number of contacts (interparticle overlaps) in the packing. Nr is the number of rattler 
particles with fewer than three contacts. The contact network is found by identifying all interparticle 
contacts, and then recursively removing rattler particles until there are none remaining in the 
packing. 

In Fig. 1 (a) we show results for the contact number zj at jamming for ellipse and dimer 
packings. We find that ellipse packings are hypostatic with zj < 2df over the range of aspect ratio 
1 < a < 2, while dimer packings are isostatic with ~ 2(ij over the same range of a. We showed 
previously [15, 20] that hypostatic packings of ellipsoidal particles possess vibrational modes that 
are quartically (not quadratically) stabilized with the number of quartic modes determined by the 
deviation from isostaticity, Ziso — zj. In contrast, all vibrational modes for dimer packings are 
quadratically stabilized since they are isostatic. 

Power-law scaling of linear shear modulus and contact number One of the hallmarks 
of the jamming transition in packings of spherical particles is the power- law scaling of the static 
linear shear modulus G and contact number z — zj with (p — (pj. Both scale as {(p — pj)^'^ for 
linear repulsive springs, which suggests that the contact number scaling controls the behavior of 
the linear shear modulus [4]. In Fig. 5, we plot G as a function oi (p — (pj for dimers (filled symbols) 
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Figure 5: Static shear modulus G versus (p — (j)j for = 480 ellipses (open symbols) and 240 dimers 
(filled symbols) at a = 1.0 (circles), 1.002 (squares), 1.01 (diamonds), 1.05 (upward triangles), 1.1 
(leftward triangles), 1.5 (downward triangles), and 2.0 (rightward triangles). The solid (dashed) 
line has slope 1 (0.5). The dot-dashed lines have the form G = 0.6((/) — <j)j)/{a — 1)*^'^^. 




Figure 6: Deviation in the contact number z from the value at jamming zj versus (f)—4>j for N = 480 
ellipses (open symbols) and 240 dimers (filled symbols) at a = 1.0 (circles), 1.002 (squares), 1.01 
(diamonds), 1.05 (upward triangles), 1.1 (leftward triangles), 1.5 (downward triangles), and 2.0 
(rightward triangles). The solid (dashed) line has slope 1 (0.5). The dot-dashed lines have the form 
z-zj = 6.3{<t>-ct>j)/{a-l)^-^'. 
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Figure 7: Shear stress Ti^y versus shear strain 7 for packings of (a) elhpses and (b) dimers at 
(j) — (f)j = 10"-*^ for several aspect ratios a = 1.0 (black), 1.05 (red), 1.1 (green), 1.2 (blue), 1.3 
(yellow), 1.4 (violet), and 1.5 (cyan). The dashed horizontal line in (b) indicates the yield stress 
for dimer packings at a = 1.3. 

and ellipses (open symbols) over a range of aspect ratios. We again find power-law scaling near 
jamming, 

G = Go{a){4>-(l)jf, (4) 

where (3 = 0.5 and Gq is weakly dependent on a for dimers. In contrast, /? = 1 for sufficiently 
small 4) — 4>j and Go(a) ~ (a — i)-0-44±o.03 £qj. giUpsgg, xhe power-law scaling is stronger for ellipse 
packings, and thus the ratio of the shear moduli Geiupse/Cdimer — > in the limit </>—)• j for all a. 
This implies that ellipse packings are much more susceptible to shear in the linear response regime. 

For jammed packings of spherical particles with linear spring interactions both G and z — zj 
scale as {(f) — (pj)^'^. We find similar behavior, G ~ z — zj, for dimer packings as a shown in Fig. 6. 
For ellipse packings, we find 

z - zj = zo{a){(j) - cl)j)^, (5) 
where zo{a) ~ (a — i)-0-35±o.i g^^^^ /3 = 1 for sufficiently small (j) — (pj. Thus, G and z — zj have the 
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Figure 8: Same stress-strain relations in Fig. 7 for (a) ellipses and (b) dimers except the shear 
stress and strain have been scaled by and S^/G, respectively. The inset to panel (a) gives the 
average value of the slope (Gint / G) of the scaled stress-strain relation for shear stresses T,xy < 
at aspect ratios a > 1.1. 

same power-law scaling with (j) — (pj even for hypostatic packings. We argued previously that the 
novel power-law scaling exponent for G and z — zj in ellipse packings originates from the quartically 
stabilized vibrational modes [15]. 

Stress-strain relations The full stress-strain behavior for ellipse and dimer packings is complex; 
it is qualitatively different for ellipses and dimers and depends nontrivially on (j) — (pj and aspect 
ratio. In Figs. 7 and 9, we show the shear stress T^xy versus strain 7 for <p — = 10~^ and 10~^. 
For ellipses at <j) — = 10^^ (Fig. 7 (a)), the shear stress is roughly linear with strain until the 
shear stress plateaus at the yield stress, = Sxy(7 — > 00), which only weakly depends on aspect 
ratio and is at least a factor of 2 smaller than that for dimers (c./. Fig. 11 (a)). We can achieve an 
approximate collapse of the stress-strain data for ellipses at (j) — = 10~^ for a > 1.1 by scaling 
the shear stress by and strain by Tiy/G as shown in Fig. 8 (a). The inset to Fig. 8 (a) shows 
that the average shear modulus defined over the wide range < Tixy < is comparable to the 
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Figure 9: Shear stress T,xy versus shear strain 7 for packings of (a) elhpses and (b) dimers at 
(j) — (pj = 10"'^ for several aspect ratios a = 1.0 (black), 1.1 (green), 1.2 (blue), 1.3 (yellow), 1.4 
(violet), and 1.5 (cyan). The insets show the same data as in the main plots, except over a smaller 
range of 7. 



linear response value, G, at small strains (c./. Fig. 5) for a > 1.1. 

The behavior of the stress-strain curves for dimers ai (j) — 4>j = 10~^ is qualitatively different 
from that for ellipses as shown in Fig. 7 (b). In particular, the approach of the shear stress to the 
yield stress plateau has significant curvature similar to the behavior found for sheared packings of 
spherical particles [25]. The scaled stress-strain curve in Fig. 8 (b) emphasizes that dimer packings 
further strain soften as the aspect ratio increases. 

Fig. 9 shows the stress-strain behavior for dimers and ellipses much closer to the jamming 
transition at (j) — = 10~^. At such small compressions, ellipse packings (panel (a)) no longer 
possess such robust, sustained linear response over the range < Tixy < Ti^. Instead, the shear 
stress is first roughly linear with slope ~ G, but then stiffens on approach to the yield stress. 
However, the most striking feature of the stress-strain curves at (p — = 10~^ is the nearly perfect 
plastic response (flat shear stress versus strain) for dimer packings (panel (b)) with a > 1.1. The 
plastic regime extends for strains from the end of the linear response regime to 7p ~ 0.1. For 
7 > 7p, the shear stress grows rapidly as it approaches the yield stress, which is similar to the 
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Figure 10: Fraction of the particle contacts at strain 7 that differ from those at 7 = for (a) 
elhpses and (b) dimers at a = 1.05 (red), 1.1 (green), 1.2 (blue), 1.3 (yellow), 1.4 (violet), and 1.5 
(cyan). 



behavior in Fig. 7 (b) at (p — cpj = W~^. 

As demonstrated in the inset to Fig. 9 (a), ellipse packings do not possess this nearly plastic 
response. In the regime 7 < 7p, the shear stress is roughly linear with a shear modulus comparable 
to G. In Fig. 10, we plot the fraction of particle contacts at strain 7 that differ from those at 
7 = 0. We find that the plastic behavior in dimer packings is accompanied by a large increase 
in the number of particle rearrangement events (changes in the contact network) over the strain 
interval < 7 < 7p. Note that the largest fraction (and rate of increase over < 7 < 7p) of particle 
rearrangements occurs for dimer packings with a = 1.5, which possess the most pronounced plastic 
response. Further work is required to elucidate the particle-scale motions that cause this plastic 
response. 

We have also calculated the nematic order parameter, S = (cos[2(0 — ^o)]); where 6q is average 
orientation of the particles, as a function of shear strain. We find that the nematic order increases 
roughly linearly with 7 up to strains of 0.2 — 0.3 in both sheared dimer and ellipse packings. Thus, 
it is possible that nematic order leads to qualitatively different effects in dimer and ellipse packings 
since dimers display plastic response, while ellipses do not. 

In Fig. 11, we show the yield stress T,^ for dimers and ellipses as a function of aspect ratio at 
(a) (/) — 4>j = 10~^, (b) 10~^, and (c) 10"'^. We find that the yield stress increases with aspect ratio, 
which acts as an effective friction coefficient [26]. However, begins to level off near a* ~ 1.4, 
which is likely related to a maximum in the nematic order near a*. In contrast to the behavior at 
small shear strains, the yield stress for dimers and ellipses becomes nearly identical near jamming 
at (/) — = 10~^. Thus, measurements of the jamming packing fraction and yield stress are 
relatively insensitive to microscale geometrical features of individual particles. 

5 Future Directions 

These studies of the structural and mechanical properties of dimer and ellipse packings raise a 
number of interesting questions that will likely spur new research activity in this area. First, there 
are many features of jammed ellipse and dimer packings that are different. Most notably, ellipse 
packings are hypostatic while dimer packings are isostatic, which gives rise to novel power-law 
scaling of the structural and mechanical properties near jamming. Thus, we have shown that 
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Figure 11: Yield stress as a function of aspect ratio a at {a) 4> — = 10 ^, (b) 10 ^, and (c) 
10^'^ for ellipses (circles) and dimers (squares). 



the macroscopic jamming behavior of anisotropic particles depends sensitively on the microscale 
geometrical features of individual particles. 

However, the dependence on aspect ratio of the jamming packing fraction and yield stress, 
which strongly affect the glass transition in thermalized systems [27], is relatively insensitive to 
whether the packings are composed of dimers or ellipses. Thus, it is not clear a priori which 
structural, mechanical, and dynamical properties are sensitive to microscale particle properties. 
Thus, we encourage reinvigorated studies of atomic, colloidal, and granular systems to determine 
under what circumstances geometrical features of individual particles play an important role in 
jamming behavior and glassy dynamics. 
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A Packing fraction 

In our numerical simulations, we consider bidisperse mixtures of dimers and ellipses in which one- 
third {Ni = 1/3) of the particles are large (with a minor axis 1.4 times that of the smaller particles, 
i.e. bi = 1.46s) and two-thirds [Ng = 2/3) of the particles are small. When calculating the packing 
fraction for rigid dimers (fused disks), we do not double count the overlapping region. Thus, we 
define the packing fraction for dimers in 2D as 



'/'dimor = 2iVs7r 

For ellipses 




^ "^(a - 1) + (a - l)^Ja{2 - a) ]. (6) 



cos 




<Aellipse = iV.vra f 1 + ^^ . (7) 
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B Contact distance 



The Pcrram and Wertheim formulation for calculating the contact distance a^j between ellipses i and 
j with orientations /ij and p,j and center-to-center direction fij involves the following minimization 
procedure [18]: 



where 



and 



j_fi 

± 1 ± x(A)/ii • fij 



an = mincrfA) = min — , :, (8) 

^ ^ /i x(A) ^ imr^J • ± ■ 



= . . . ) [ ' „^ , (10) 




m = -Tz — r-:v • (11) 



Determining Amin that minimizes a{X) (Eq. 8) involves solving for the roots of a quartic polynomial 
in A for 2D bidisperse systems [20]. The polynomials can be expressed analytically in terms of /tj, 
fij, rij, and the major and minor axes of particles i and j, and then solved using Newton's method. 



C Packing-generation algorithm 

In Sec. 3, we outlined our procedure to generate static packings of dimers and ellipses. Here, we 
provide some of the numerical parameters involved in the simulations. For the energy minimization, 
we employ the conjugate gradient technique [23], where the particles are treated as massless. The 
two stopping criteria for the energy minimization are Vt — Vt-i < Hoi = lO"-*^^ and Vt < Vmm = 
10""^^, where Vt is the potential energy per particle at iteration t, and the target potential energy 
per particle of a static packing is Vtol < V/N < 2Vto\- For the first compression or decompression 
step we use the packing- fraction increment = 10~^. Each time the procedure switches from 
expansion to contraction or vice versa, Acj) is reduced by a factor of 2. Using the packing generation 
procedure with these parameters, we are able to locate the jamming threshold in packing fraction 
to within 10~^ for each static packing. 



D Calculation of forces and torques 

In this appendix, we provide specific details for calculating the interparticle forces and torques for 
dimers and ellipses, which arc required to perform energy minimization and evaluate the shear 
stress. The forces and torques can be obtained from the interaction potential (Eq. 1) using gener- 
alizations of Fij = dV/dfij, where F^j is the force on particle i due to particle j. 
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Figure 12: Definition of the point of contact Cij for ellipses i and j that are (a) 'just touching' and 
(b) overlapped, pij is location of the point of contact relative to the center of mass of ellipse i. In 
the overlapped case, the effective point of contact C^j is given by p^^- = Pij{rij / Uij). 



Dimers For dimers, the interaction force on monomer ki belonging to dimer i from monomer kj 
belonging to a distinct dimer j is 

F =J^_ (12) 

' dfkk 

The total force on dimer i is obtained by summing over all monomers hi belonging to dimer i, all 
dimers j different from i, and all monomers kj belonging to dimer j: 

hi 'J kj 

The torque on dimer i arising from an interaction between monomer ki on dimer i and monomer 
kj belonging to dimer j is given by 

Tk„k, =rki X Fk,,kj, (14) 

where r^. = di{cos6iX — sin 0iy) is the vector from the center of dimer i to the center of monomer 
ki, di = bi{a — 1) /2, and 9i gives the orientation of dimer i. The total torque on dimer i, Ti, is 
obtained by summing Tk^^kj over all monomers ki on dimer i, all dimers j distinct from i, and all 
monomers kj on dimer j. 

Ellipses For ellipses i and j, the interparticle force depends explicitly on how the contact distance 
(Tjj varies with the vector separation fij-. 

dV f Cij Tij dcFij 



dvij \rij Gij diijj 

where ^ = x,y. To calculate the torque, one must specify the point of contact. For 'just touching' 
ellipses i and j are in contact at only one point, as shown in Fig. 12 (a), the location pij of the 
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point of contact Cij (relative to the center of mass of ellipse i) is unambiguous and given by 

Pij = Pij{cos{ipij + ei)x + sm{'iljij + ei)y) (16) 

+ 



ai 



iaxiihij = jr^> (18) 

' «'l + ^?tan(e,-^.)||- 

where cos Pij = fii ■ rij, cos @ij = xi ■ rij. The torque on ellipse i from j is then 

As shown in Fig. 12 (b), upon compression, ellipses are no longer 'just touching', and thus Eq. (16) 
for the point of contact Cij is no longer exact. In this case, we scale pij by Vij/aij, which yields an 
effective point of contact C^j that is within the overlap region of the two ellipses. 



E Calculation of shear stress 

For systems composed of spherical particles, the correct form for the stress tensor in 2D, where 
a, P = x,y, is the virial expression [24]: 

1 ^ 

= ^ Yl (Pijo'rijp + Fijfsnja) , (20) 
i>j=l 

where Fija is the a-component of the force Fij on particle i arising from an overlap with particle j, 
Tijp is the /3-component of the vector fij from the center of mass of particle j to that of particle i. 

The correct form for the stress tensor in 2D for systems composed of anisotropic particles 
is the Love expression [16]: 

1 ^ 

= ^ {FijaPijIi + FijfiPija) , (21) 

where pjj^ is the /3-component of the vector from the center of mass of particle i to the point 
of contact Cij with particle j. Note that the Love expression reduces to the virial expression for 
spherical particles. 

In our studies of simple shear, we focus on the off-diagonal component of the stress tensor 
Tjxy Calculating the point of contact at each shear strain is computationally expensive; we have 
therefore used the virial expression instead of the Love expression to quantify the shear 
stress for both dimer and ellipse packings. As a check, we measured both "E^y and "E^y for dimers 
as a function of aspect ratio and compression. Fig. 13 shows that they give quantitatively similar 
results for a = 1.1 and a = 1.5 (for 7 < 0.2), and qualitatively similar results for a = 1.5 at large 
strain. In particular, the plastic response of dimer packings at small compressions is unaffected by 
the choice of the definition of the shear stress. 



14 




References 

[1] A. J. Liu and S. R. Nagel, "Jamming is not just cool anymore," Nature 396, 21 (1998). 

[2] Z. Zhang, N. Xu, D. N. Chen, P. Yunker, A. M. Alsayed, K. B. Aptowicz, P. Habdas, A. J. Liu, 
S. R. Nagel, and A. G. Yodh, "Thermal vestige of the zero-temperature jamming transition," 
Nature, 459, 230 (2009). 

[3] M. van Hecke, "Jamming of soft particles: geometry, mechanics, scaling and isostaticity," J. 
Phys.: Condens. Matter 22, 033101 (2010). 

[4] C. O. O'Hern, L. E. Silbert, A. J, Liu, and S. R. Nagel, "Jamming at Zero Temperature and 
Zero Applied Stress: The Epitome of Disorder", Phys. Rev. Lett. 68, 011306 (2003). 

[5] W. G. Ellenbroek, E. Somfai, M. van Hecke, W. van Saarloos, "Critical Scaling in Linear 
Response of Frictionless Granular Packings near Jamming", Phys. Rev. Lett. 97, 258001 (2006). 

[6] L. E. Silbert and A. J. Liu and S. R. Nagel, "Vibrations and diverging length scales near the 
unjamming transition," Phys. Rev. Lett. 95, 098301 (2005). 

[7] W. G. Ellenbroek, M. van Hecke, and W. van Saarloos, "Jammed frictionless discs: Connecting 
local and global response," Physical Review E 80, 061307 (2009). 

[8] M. Wyart and S. R. Nagel and T. A. Witten, "Geometric origin of excess low-frequency 
vibrational modes in weakly connected amorphous solids," Europhys. Lett. 72 486 (2005). 

[9] A. V. Tkachenko and T. A. Witten, "Stress propagation through frictionless granular materi- 
als," Phys. Rev. E 60, 687 (1999). 



15 



[10] K. Desmond and S. V. Franklin, "Jamming of 3D prolate granular materials," Physical Review 
E7S, 031306 (2006). 

[11] N. Mori and R. Semura and K. Nakamura, "Simple shear flows of suspensions of Brownian 
ellipsoids interacting via the Gay-Berne potential," Molecular Crystals and Liquid Crystals 
367, 445 (2001). 

[12] A. Donev, R. Connelly, F. H. Stillinger, and S. Torquato, "Underconstrained jammed packings 
of hard ellipsoids," Phys. Rev. E 75, 051304 (2007). 

[13] Z. Zeravcic, N. Xu, A. Liu, S. Nagel, W. van Saarloos, "Excitations of ellipsoid packings near 
jamming," Europhys. Lett. 87, 26001 (2009). 

[14] A. Donev, I. Cisse, D. Sachs, E. A. Variano, F. H. Stillinger, R. Connelly, S. Torquato, and P. 
M. Chaikin, "Improving the density of jammed disordered packings using ellipsoids," Science 
303, 990 (2004). 

[15] M. Mailman, C. F. Schreck, C. O. O'Hern, and B. Chakraborty, "Jamming in systems com- 
posed of frictionless ellipse-shaped particles", Phys. Rev. Lett. 102, 255501 (2009). 

[16] S. F. Edwards and D. V. Grinev, "Transmission of stress in granular materials as a problem 
of statistical mechanics," Physica A 302, 162 (2001). 

[17] D. J. Cleaver, C. M. Care, M. P. Allen, M. P Neal, "Extension and generalization of the 
Gay-Berne potential", Phys. Rev. E 54, 559 (1996). 

[18] J. W. Perram and M. S. Wertheim, "Statistical mechanics of hard ellipsoids I. Overlap algo- 
rithm and the contact function", J. Comp. Phys. 58, 409 (1985). 

[19] J. W. Perram, J. Rasmussen, E. Praetgaard, "Ellipsoid contact potential: Theory and relation 
to overlap potentials," Phys. Rev. E. 54, 6565 (1996). 

[20] C. F. Schreck and C. S. O'Hern, "Computational methods to study jammed systems," in 
Experimental and computational techniques in soft condensed matter physics. Edited by J. S. 
Olafsen (Cambridge University Press, London, 2010). 

[21] G.-J. Gao, J. Blawzdziewicz, and C. S. O'Hern, "Studies of the frequency distribution of 
mechanically stable disk packings", Phys. Rev. E 74, 061304 (2006). 

[22] G.-J. Gao and J. Blawzdziewicz and C. S. O'Hern, "Geometrical families of mechanically stable 
granular packings," Phys. Rev. E 80, 061303 (2009). 

[23] W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes in 
Fortran 77 (Cambridge University Press, New York, 1986). 

[24] M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, 
New York, 1987). 

[25] N. Xu and C. S. O'Hern, "Measurements of the yield stress in frictionless granular systems," 
Phys. Rev. E 73, 061303 (2006). 

[26] L. E. Silbert, D. Ertas, G. S. Crest, T. C. Halsey, D. Levine, and S. J. Plimpton, "Granular 
flow down an inclined plane: Bagnold scaling and rheology," Phys. Rev. E 64, 051302 (2001). 



16 



[27] C. De Michclc and R. Schilling and F. Sciortino, "Dynamics of uniaxial hard ellipsoids," Phys. 
Rev. Lett. 98, 265702 (2007). 



17 



